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Primordial non-Gaussianity is a crucial test of inflationary cosmology. We consider the impact 
of non-Gaussianity on the ionization power spectrum from 21 cm emission at the epoch of reion- 
ization. We focus on the power spectrum on large scales at redshifts of 7 to 8 and explore the 
expected constraint on the local non-Gaussianity parameter /nl for current and next-generation 21 
cm experiments. We show that experiments such as SKA and MWA could measure /nl values of 
order 10. This can be improved by an order of magnitude with a fast-Fourier transform telescope 
like Omniscope. 
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Introduction. An inflationary epoch in the early 
universe [J, Q has been established as a solution to 
the cosmological horizon and flatness problems over 
the past three decades, most recently through high- 
precision measurements of the cosmic microwave back- 
ground (CMB) by the Wilkinson Microwave Anisotropy 
Probe (WMAP) [3|. The inflationary hypothesis predicts 
an epoch of exponential growth lasting at least 60 e-folds 
resulting in almost Gaussian scale-invariant density per- 
turbations 3- 

A powerful mechanism to distinguish between infla- 
tion models is the amplitude and scale dependence of 
mild non-Gaussianity in perturbations of the primordial 
density field. Canonical single field inflation models pre- 
dict primordial non-Gaussianity (bispectrum) of the local 
form |/nl| ^ 1 [1, 0|, while evolution after inflation gen- 
erates non-local bispectrum with effective /nl = 

The best current constraints of ±25 on local 
/nl [13, El are from WMAP data. A future measure- 
ment of /nl = C'(l) could reveal the existence of physics 
beyond the standard single field slow-roll inflationary sce- 
nario. 



quantify the influence of non-Gaussianity of the local 
form on the 21 cm power spectrum, and then test this 
via numerical simulations of the ionization distribution. 



We show that radio interferometric probes jl2| - |l5l | of 
21 cm emission from spin-flip transitions of neutral hy- 
drogen at the epoch of reionization (EoR) [3l can result 
in constraints on /nl at the same level as Planck IVf, 
and less than unity in the most optimistic experimental 
proposal. Previous studies have explored primordial non- 
Gaussianity in the bispectrum of ideal 21 cm experiments 
prior to the EoR 18|, ll9| . In this work, we consider scale 
dependent bias in the power spectrum of ionized hydro- 
gen resultiirg from departures from Gaussian initial con- 
ditions 2^ 21 1 . Our constraints from 21 cm emission do 
not require an ionization-clean cosmology, i.e., a priori 
knowledge of the spectrum of fluctuations in the ionized 
fraction. 

The rest of the letter is arranged as follows. We flrst 



We review the assumed noise properties of LOFAR [l^ , 



MWA [IJI, SKA [13[, and Omniscope [151, and forecast 
constraints on /nl based on a Fisher matrix analysis. For 
these forecasts, we fix the parameters of our fiducial flat 
ACDM model to agree with WMAP7 

Effect of Non-Gaussianity on the 21 cm Power 
Spectrum. We decompose the 21 cm power spectrum 
at redshift z in terms of its angular dependence 23| , given 
by /i = k • fi = cos(6'), where 9 is the angle between 
wavevector k and line of sight (LOS) vector n: 

PAT(k, z) = 3^ss{k, z) ~ 2^,s{k, z) -f ^,,(/c, z) 
+2 [^55(fc, z) - ^,5(fc, z)] fi^ + ^ssik, z)^i\ (1) 

We deflne ^ss = f'^x^Pss [3, where Pss is the lin- 
ear matter power spectrum, numerically obtained from 
a modified version of CAMB [2^, is the mean neu- 
tral fraction of hydrogen such that the ionized fraction 
X{ = 1 — xh, and Tf,(z = 7.5) ~ 0.026 K is the spatially 
averaged brightness temperature. We consider only large 
enough scales {k < 0.15/Mpc) such that the ionization 
power spectrum J3^xx — ^ss and the ionization-density 
cross spectrum — bx^ss, where is the bias of ion- 
ized regions. Our numerical simulations in Fig. [1] show 
that this is an excellent approximation. 

We define u as the Fourier dual of = diCi + Ojij + 
A/gfc, where 9i and 9j encode the angular location on 
the 2D sky, and A/ measures the difference in frequency. 
The 21 cm power spectrum is extended to u-spacc in 
which measurements arc made: 



PAr(u,z) = PAr(k,2:)/(x'(z)2/(z)), 



(2) 



where xi^) is the comoving distance to a given redshift, 
y{z) = A2i(l + zY /H{z) translates between intervals in 
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FIG. 1. Ionization power spectra with non-Gaussianity 
of the local form from numerical simulations. We show 
/nl = (0, 20, 100) (dot-dashed, dashed, solid) for efficiency 
^ — (5.8, 3.0) (thin black, thick red) at 2 = 7.5, where 
xu = (0.50,0.75). For /nl = 100 cases, sample variance from 
simulations is in form of green bands about the mean, and 
analytical fits corresponding to = 1 are in dotted lines. 



frequency and distance, and A21 = A(z)/(l + z) = 0.21 ni. 
We convert between u and k spaces via u±_ — x{z)k± = 
2ttL/X{z), where L is the baseline, and U|| = y{z)k^\. 

Given non-Gaussianity of the local form, Bardeen's 
gauge invariant potential field $ is related to a pure gaus- 
sian random field (j) at nonlinear order 0, [1^ : 



$ng(x)=0(x) + /nl (0'(x)-((/.2)) 



(3) 



In the high-peaks formalism /nl infiuences biased tracers 
of the underlying matter distribution as a scale depen- 
dent correction to the large scale bias [l^llll . This enters 
as ^^s/^5S = b^ + Ab^, ^:r. J^S5 = {bx + A6^)^ with 

A6,(fc, z) = 3(6, - l)hi^n^Hl5Bl{D{z)eT{k)), (4) 

where Hq is the Hubble constant, $7^ is the present den- 
sity parameter of matter, D{z) is the linear growth func- 
tion of density perturbations, and T(fc) is the transfer 
function relating present and primordial power spectra. 
The quantity 5b is the average critical collapse density of 
HII regions [27|. We leave the bias as a free parameter, 
although hx, 5b, and xh would all be related in a given 
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FIG. 2. Number density of baselines for LOFAR (solid), SKA 
(dotted), MWA (dashed), and Omniscope (dot-dashed). 

^periment A^ant i^min (m) FOV (deg^) Ae(m^) 
LOFAR 32 
MWA 500 
SKA 1400 
Omniscope 10® 

TABLE I. Experimental specifications for the telescopes. The 
antenna number only account for those inside the nucleus and 
core (e.g. for SKA we use 1400 of a total 7000 antennae). 
The system temperature Tsys = 390 K, bandwidth is 6 MHz, 
observation time is 4000 hours, and effective area at z = 7.5. 



model of reionization. The scale-dependence of the bias 
in Abx is clearly evident in the ionization spectra from 
our simulations in Fig. [1] We find that 5_b ~ 1 fits the 
large-scale /nl induced rise to the ionization spectrum. 

Numerical Simulations with Non-Gaussian Ini- 
tial Conditions. We perform simulations of the ion- 
ization distribution during the EoR for /nl = (0, 20, 100) 
and ionization efficiency C = (3.0,5.8), in a box of co- 
moving length 3000 Mpc, with a modified version of 
SimFast21 [2^. The initial matter density field is 
computed from the Poisson equation with non-Gaussian 
gravitational potential $NG(k). We show the spectra 
from these simulations in FigfU from which 5^ ^ 1. 

We compare this result to the theoretical prediction. 
The critical density for collapse of an ionized region of 
mass m is obtained from the collapse fraction /coii (27j : 



5B{m,z) = 5c - K{C)\/2 {a'^{mnun,z) -a^{m,z)), (5) 

where 5c ~ 1.68 is the critical collapse density of mat- 
ter, (T^(m, z) is the variance of the density fluctuations, 
and mmin corresponds to a virial temperature of 10^ K. 
Moreover, K{() = erf~^(l — C~"^), where C = 'Tiion/'Tigai 
is the ionization efficiency [27|. We evaluate 5b as an 
average over the fraction of space filled by HII bubbles 



as in Ref. [27[. Given this prescription, we find 5b — 1.1 
(less than (5c as C > 1), matching the simulation results 
well. This becomes Sb = 1.2 if we only average over the 
mass function. For simplicity, we fix (5b = 1. 

As noted earlier, the bias bx, collapse threshold 5b, and 
xh are expected to be interrelated in a given reionization 
scenario. This is evident in Fig. [1] where we see that a 
factor of 2 change in ( changes the bias by about 15%. 
This change is subdominant to the impact of xjj (hnear 
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function of C) on the 21 cm power spectrum. In a more 
optimistic scenario, one could envision constraining xn 
(or Q together with /nl without as a free parameter. 
We also considered the impact of variations in xh and 
/nl on 5b- Changing xn by a factor of two only affects 
6b by 8% given (1 - xh) = C/coii- Nonzero /nl skews 
5B{m, z) through its influence on /cou. Using the results 
of Ref. [30|, we estimate is only perturbed by 4% even 
for /nl = 100. This is because the sensitivity to /nl in- 
creases with mass, while the mass scales that contribute a 
majority of the integral over the mass function lie within 
an order of magnitude of the minimum halo mass. 

21 cm Noise Power Spectrum. The noise power 
spectrum of 21 cm fluctuations is expressed as p3. [sH 

PN{uuz)^{\^{z)T,y,{z)/Ae{z)f Hton{u^)), (6) 

where the sky-dominated system temperature Tgys — 
280((l + z)/7.4)^-^K [ail, h is the total observation 
time, and Ae{z) cx A^(z) is the effective collecting area 
(listed in Table|I|. Here, n(u±) encodes the number den- 
sity of baselines shown in Fig. [2l computed as the auto- 
correlation of the array density for each of the surveys. 

The array distributions are composed of a nucleus with 
full coverage fraction and a core with power law r~^. The 
nucleus radius is i?„ = ■\/ ?7^ant/(po7r), where po is the 
2D array density of the nucleus, and A^ant is the number 
of antennae of each experiment (see Table |I| . The core 
radius is by construction Rc = Rn exp((l — r])/{2r])) [3]. 
The most optimal choice of ij for constraints on /nl de- 
pends on the particular experiment and bandwidth B 
considered, but for comparison with prospective con- 
straints on other cosmological parameters in Table V of 
Ref. (ill, we choose rj = 0.8 for [LOFAR, MWA, SKA] , 
whereas all of Omniscope's antennae lie in the nucleus. 

We assume residual foregrounds can be ignored beyond 



k\\ > 2Tr /{yB) [3l|, but also consider the case where fore- 



grounds can be removed on larger scales (Fig [3]). 

Fisher Matrix Forecasts. We evaluate the prospec- 
tive constraints on /nl from the 21 cm power spectrum 
at the EoR via the Fisher matrix formalism. The sum- 
mation involves pixels in (k±,k\\) of thickness (e^,e||) ~ 
(Afcx/fci,Afc,|/fc||) = (0.1,0.1): 
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(7) 

We have verified that our forecasts are robust to vari- 
ations in the step sizes of parameter space and k-space. 
The measurement error consists of the sum of the sarnple 
variance and thermal detector noise over half-space [31| : 



SPat{u) = (Pat(u) + PN{u^))/VN^r 



(8) 



The number of modes falling in each pixel is given by 
Nm = 27rfcj_Afcj^Afc||V(z)/(27r)^, such that the volume 
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FIG. 3. Top: Marginalized /nl constraints for cases with 
noise (thick) and without noise (thin), which overlap for Om- 
niscope. We consider a bandwidth of 6 MHz, but assume fore- 
grounds can be removed on scales larger than fcy = 2n /{yB). 
Bottom: Marginalized /nl constraints as function of band- 
width and number of antennae. The bandwidth limits the 
number of modes and largest scale probed along the LOS 



(via the survey volume V x B and k]l 



1/B), whereas 



a larger number of antennae for fixed array density increases 
the survey resolution and number of perpendicular modes (via 
n{u±), on large scales oc A''aiit, and it™'"' oc y/N^It). The color 
coding is the same as for the top panel. 



sampled V{z) = x^yB x FOV, where FOV denotes the 
field of view of the telescope (often equal to /A^). 
For a single redshift bin at z = 7.5, we fiducially let 
= 2.3 and x-a = 0.5. The bandwidth B = 6 MHz 
limits fcfp'" = 2TT/{yB) > 0.063/Mpc [3l|, and non- 
linearities force k]^^^ ^ 2/Mpc. The ranges in k± at 
the central redshift are [0.039, 0.25] /Mpc for LOFAR, 
[0.0016, 0.040] /Mpc for MWA, [0.0039, 0.17] /Mpc for 
SKA, and [3.9 x 10""*, 0.44] /Mpc for Omniscope. How- 
ever, due to our narrow focus on /nl at the largest scales 
in which the 1/fc^ boost becomes significant, in practice, 
we only keep modes up to fc™^'' = 0.15/Mpc. 



Results. In quantifying our constraints on /nl, we 
fix the underlying cosmology. By only considering large 
enough scales for which the ratio of the ionization and 
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matter spectra is constant in a universe without non- 
Gaussianity, the free parameters in a single redshift bin 
at in = 0.5 are hmited to (/nl, ^k, ^h)- With Planck 
priors on the standard cosmological parameters [l7j . in 
particular the matter power spectrum normalization A|j, 
cold dark matter density Oc/i^, spectral index Ug, and 
its running dnj/dln/c, we find the /nl constraints from 
[LOFAR, MWA, SKA] are robust to the assumption of a 
fixed cosmology at the 10% level, while the same level of 
robustness for Omniscope is achieved after including its 
constraints on (ris, dris/d In fc, fic/i^) from small scales. 
The constraints on /nl will depend on the fiducial b^, 
but we do not explore this issue here. 

Fig. |3] (top) shows /nl constrained as function of the 
minimum LOS wavenumbcr, limited by the experimen- 
tal ability to remove foregrounds. Imposing fcj"'" = 
27r/(y_B) = 0.063/Mpc [3l|, we find the constraints for 
[LOFAR, MWA, SKA, Omniscope] are equal to cr(/NL) = 
[700, 100, 50, 4], which reduces to cr(/NL) = [100, 30, 40, 4] 
when instrumental noise is neglected. These constraints 
improve for telescopes with increased ability to probe 
larger LOS scales. When arbitrarily large scales along 
the LOS can be probed, we find cr(/NL) = [200, 6, 10, 0.6], 
which reduces to ^(/nl) = [70,5,10,0.6] when noise is 
neglected. The constraints plateau for fc]"'" due to 
the nonzero fc™™ set by the minimum experiment base- 
line. As fc™™ decreases, our assumed MWA configuration 
becomes somewhat better than the SKA configuration in 
constraining /nl due to its smaller minimum baseline, 
allowing larger scales to be probed by the telescope. 

In Fig. [3] (bottom), we consider a minimum LOS scale 
set by fcJi"'" — 2TT/{yB), but allow an order of magni- 
tude variation in bandwidth and telescope antenna num- 
ber. The bandwidth is inversely proportional to the mini- 
mum LOS wavenumbcr and linearly increases the volume 
probed, whereas larger number of antennae for fixed ar- 
ray density increases the maximum baseline as y/ A^ant 
and linearly boosts the baseline density (thereby decreas- 
ing the noise). The contours show increased bandwidth 
is more powerful in the search for /nl, in particular for 
SKA and Omniscope that have small instrumental noise. 
This is because their signal-to-noise is already close to 
the cosmic variance limit, and our power spectrum cutoff 
at fc = 0.15/Mpc makes us insensitive to the increasing 
number of small scale modes. Extending the considered 
modes to scales of fc = 2/Mpc (incorporating modeling 
of the exponential tail with very strong priors on the new 
free parameters) improves the constraints by up to factor 
of 2 for the different experimental configurations. 

We have also considered the case where the bias and 
ionization fraction are fixed. In this scenario, the /nl 
constraints improve by a factor of 1.5 up to a factor of 10 
for the various cases and experiments considered. For the 
fiducial configurations alone, the /nl constraints improve 
by factors of 2 (MWA) to 3 (LOFAR, SKA, Omniscope) 



when fixing the bias to be a function of the ionization 
fraction. When only information from scales larger than 
^max ^ 0.10/Mpc is available (compared to 0.15/Mpc 
assumed throughout the paper), the constraint on /nl 
degrades by up to a factor of 2 when marginalizing over 
bx and xh, and by up to 30% when bx and xh are fixed. 

Conclusions. The search for a signature of primor- 
dial non-Gaussianity is a key test of inflationary theo- 
ries. Large values for the non-Gaussianity parameter, 
/nl ^ 1, will rule out standard single field inflationary 
models. We have considered the impact of primordial 
non-Gaussianity on the ionization power spectrum from 
21 cm emission at the epoch of reionization, which pro- 
vides an alternative approach to constrain /nl relative to 
the cosmic microwave background and large-scale struc- 
ture. Wc flnd that /nl can be constrained to an accuracy 
of order 10 with future 21 cm telescopes like SKA and 
MWA. This improves by an order of magnitude for a 
fast-Fourier transform telescope like Omniscope, thereby 
opening a new window to inflationary physics. 
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